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We have used time-of-flight inelastic neutron scattering to measure the spin wave spectrum of 
the canonical half-doped manganite Pro.sCao.sMnOa in its magnetic and orbitally ordered phase. 

The data, which cover multiple Brillouin zones and the entire energy range of the excitations, are 
compared with several different models that are all consistent with the CE-type magnetic order, 
but arise through different exchange coupling schemes. The Goodenough model, i.e. an ordered 
state comprising strong nearest neighbor ferromagnetic interactions along zig-zag chains with anti¬ 
ferromagnetic inter-chain coupling, provides the best description of the data, provided that further 
neighbor interactions along the chains are included. We are able to rule out a coupling scheme in¬ 
volving formation of strongly bound ferromagnetic dimers, i.e. Zener polarons, on the basis of gross 
features of the observed spin wave spectrum. A model with weaker dimerization reproduces the 
observed dispersion but can be ruled out on the basis of discrepancies between the calculated and 
observed structure factors at certain positions in reciprocal space. Adding further neighbor interac¬ 
tions results in almost no dimerization, i.e. recovery of the Goodenough model. These results are 
consistent with theoretical analysis of the degenerate double exchange model for half-doping, and 
provide a recipe for how to interpret future measurements away from half-doping, where degenerate 
double exchange models predict more complex ground states. 


I. INTRODUCTION 

A characteristic of many materials exhibiting strong 
electronic correlations is a delicate balance between mul¬ 
tiple competing phases which can often be tuned with 
relatively modest changes in external parameters such 
as temperature, magnetic field, pressure or chemical 
doping^. Preeminent amongst strongly correlated elec¬ 
tron systems are transition metal oxides, which dis¬ 
play phenomena as diverse as superconductivity^, orbital 
order—, and complex magnetic order—. Most families of 
transition metal oxides have been studied extensively, 
and manganites are no exception. 

Interest in the manganites in recent times has broadly 
centered on two phenomena - colossal magnetoresistance 
(CMR)^ and multiferroicity^. CMR manganites are in 
many ways typical of strongly correlated electron sys¬ 
tems, in that they exist at the cusp of several compet¬ 
ing phases^. For the particular case of (nearly) cubic 
Ri_a;Aa;Mn03 (where R is a rare-earth or La, and A is Ca, 
Sr, Ba, or Pb) and its layered analogs (R,A)2Mn04 and 
(R,A)3Mn207 such phases include charge-ordered, an¬ 
tiferromagnetic, ferromagnetic metal, and ferromagnetic 
insulator. In these perovskite manganites each Mn ion 
has three electrons in the t 2 g d-orbitals, and a fraction 
(1 — x) have one in the higher energy Cg orbitals, whose 
spins are aligned parallel to the local t 2 g moment by 


strong intrasite exchange. The essential physics is a com¬ 
petition between delocalization of the eg electrons, favor¬ 
ing a homogeneous, fully spin polarized, ferromagnetic 
metal via double exchange, and localization due to lattice 
distortions trapping the eg electrons to form polarons^. 
For X ~ 0.3 the materials are typically ferromagnetic 
metals, with a sea of polarons in the insulating paramag¬ 
netic phase. The most marked CMR effects often occur 
for those materials in which the polarons order in the half 
doped state, x ~ 0.5, to give a charge and orbitally or¬ 
dered insulating state for T < Tqo, the charge-ordering 
temperature. In the case of Pro.5Cao,5Mn03(PCMO) 
Tco = 260 K and the transition from a paramagnetic 
to antiferromagnetic phase occurs at Tn = 180 K^. 

Despite its relevance to the defining properties of the 
CMR manganites, the nature of the charge-ordered state 
is not yet well-understood. At half-doping, the pic¬ 
ture widely assumed is of CE-type charge and mag¬ 
netic ordering, illustrated in Fig. UKa), first proposed 
by Goodenougbi^. In this original picture, half the Mn 
ions are Mn®+ (hence spin 2) and half are Mn"*^+ (spin 
3/2) in a plane of the pseudocubic lattice; the coupling 
is ferromagnetic (FM) between ions along the zig-zag 
chains, with antiferromagnetic (AFM) inter-chain cou¬ 
pling. Along the zig-zag chains the nearest neighbor Mn 
sites are occupied by alternating Mn®’*' and Mn"^’*' ions, 
with the lobes of the d3^2_,,2-type orbitals on the Mn®+ 
sites pointing at the Mn^+ ions on the zig-zags. These 
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planes are then stacked antiferromagnetically. There is 
considerable experimental evidence in favor of the Good- 
enough model. Nnmerous diffraction and x-ray resonant 
scattering results point to two inequivalent Mn sites in 
various half-doped manganites^i^i^— , and the CE-type 
magnetic structure is universally accepted^ d^d^ . The 
charge modulation is markedly less than expected in this 
simple picture in various half-doped manganites, how¬ 
ever, including Lao,5Sri,5Mn04, Ndo.sSi'o.sMnOa and 
Prg.sCao.sMnOa. Bond valence sums applied to the Mn 
ions typically reveal Mn valences of approximately 3 . 4 - 1 - 
and 3 . 6 - 1 - rather than the expected 3 - 1 - and 4 + at the 
two sites^di, and X-ray resonant scattering at the Mn 
K-edge also suggests that the valences differ by the same 

amoun t^^d^ . 

The formation of the charge-ordered and CE-type mag¬ 
netic structure at half-doping can be understood theoreti¬ 
cally in terms of the degenerate double exchange (DDEX) 
mode b^d^ . The DDEX model is a generalization of dou¬ 
ble exchange that allows for different occupancies of the 
two-fold degenerate Cg orbitals on a lattice of Mn sites, 
finite Hund’s coupling between the eg and t2g orbitals, 
and superexchange between the core moments of the Mn 
sites. Frustration between double exchange, which fa¬ 
vors delocalization and ferromagnetic alignment of the 
Mn core moments, and the antiferromagnetic superex¬ 
change results in a cooperative orbital and charge order¬ 
ing that favors ferromagnetism along zig-zag chains and 
antiferromagnetic stacking of the chains. 

At half-doping and for suitable choice of (renormal¬ 
ized) superexchange and electron hopping energies, the 
DDEX model predicts site-centered charge ordering, with 
charge disproportionation of ^ O.le which is similar to 
the experimental values. Below half-doping there is a 
range of hole doping, x, over which ferromagnetic dimer 
formation takes place in the ground state. The limiting 
case of bond-centered charge ordering, in which the Mn 
ions in a dimer are equivalent, occurs for x ^ 0.4 . This 
limiting case of a dimer resembles the Zener polaro n^^i^° 
(ZP), in which adjacent Mn ions pair up via delocaliza¬ 
tion of the lone Cg electron to produce a single strongly 
bound unit with spin 7 / 2 . The Zener polarons form a 
herringbone pattern, as shown in Fig. Hb), with AFM 
coupling between parallel units, and EM coupling be¬ 
tween perpendicular units. 

The idea of the ZP came to prominence following a 
key single crystal neutron diffraction experimental on 
Pro,6Cao.4Mn03 performed at a temperature below Tco 
in which the orientation of the elongated MnOe octahe- 
dra and the two similar (although not identical) Mn sites 
in the structural refinement suggest the formation of Mn- 
Mn pairs. The concept of Zener polarons also received 
theoretical support from ab initio calculations performed 
at the time in half-doped manganites^^— . However, 
the existence of Zener polarons has proved controversial. 
Notwithstanding the evidence in support of the Goode- 
nough model for half-doped manganites cited above, X- 
ray resonant scattering results point to two inequivalent 


Mn sites in Pro,6Cao.4Mn03 itself, which in the full 
analysis is shown not to be consistent with the form of 
ZP defined in ref. [U Indeed, of the studies cited in favor 
of the Goodenough modelll^— , many explicitly rule out 
the existence of Zener polarons. On the other hand, re¬ 
ports of high resolution transmission electron microscopy 
and electron diffraction experiments claim to confirm 
the existence of ZP-type ordering in Pri_a,Gaa;Mn03, 
X = 0.3 — 0 . 5 , with phase coexistence of the ZP and 
the conventional charge-ordered structure2Si21. 

The existence of two inequivalent sites discussed above 
as evidence against ZP formation in itself does not pre¬ 
clude the existence of dimers made from adjacent Mn 
ions. There are a continuum of possibilities for hole dop¬ 
ing between the ZP (bond centered) and the Goodenough 
(site centered) limits, in which the two Mn sites in the 
ferromagnetic dimers that the DDEX model predicts are 
inequivalent, with different occupancy of the two eg or¬ 
bitals at the two positions. The dimers break inversion 
symmetry and this leads to the possibility of ferroelectric- 
ity in manganites; indeed there have been reports^^*^ of 
spontaneous electric polarization in Pri_3;Gaa;Mn03 for 
a range of x, including x = 0 . 5 . Motivated by the possi¬ 
bility of multiferroicity, more sophisticated theories that 
extend the DDEX model to incorporate on-site Coulomb 
repulsion together with MnOg tilting, and / or Jahn- 
Teller distortions have been developed, which have been 
analysed in terms of model calculations or a combination 
of these with DFT for half-doped manganites^SiSi. These 
approaches reveal instability towards some degree of Mn- 
Mn dimerization for a range of parameters in the model, 
as do full ab initio calculations^^*^ for Pro,5Cao.5Mn03. 

Because the origin of the ferromagnetic interaction is 
double exchange, the characteristic strength of the mag¬ 
netic exchange can reasonably be expected to be that 
typical of ferromagnetic manganites^^— , which in turn 
is consistent with estimates made from realistic values 
for the band width o^^'^^ . The values of the magnetic ex¬ 
change constants within and between dimers will in gen¬ 
eral be different because of the inequivalent eg orbital 
occupancies on the two Mn sites. In the extensions to 
the DDEX model that incorporate interaction with the 
lattice, Jahn-Teller distortions (which couple with the Cg 
orbital occupancies) and alternating modulation of Mn- 
0 -Mn bond angles also lead to different intra and inter 
dimer exchange. 

The motivation for the current work is two-fold. First 
is the question of whether or not the Mn ions are paired 
in a dimer-like fashion in Pro.sGap.sMnOs. On the 
basis of previous inelastic neutron scattering measure¬ 
ments of other half-doped manganite s^^d^d^ we might 
expect the Goodenough model to provide a better de¬ 
scription than models involving dimerization. However, 
since Pri_2,Gaa;Mn03 is the material in which the possi¬ 
ble existence of Zener polarons came to prominence, we 
believe that it is important to confirm whether or not 
this is the case. This then leads to our second motiva¬ 
tion, which is to present an exhaustive account of how 
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FIG. 1 . (Color online) (a) Illustration of the Goodenough model (see text). The spins in different one-dimensional chains 
(zig-zags) are indicated by either red (dark) or purple (light) arrows. The orbitals on the Mn 3 -|- sites are d^^2_^2 or d^y2_^2 
in character, where x and y point along the two perpendicular Mn-Mn directions along the zig-zags. In the model as orig¬ 
inally proposed by Goodenougbi^ there is a ferromagnetic interaction between nearest neighbors along a zig-zag, Jfi, and 
antiferromagnetic coupling between nearest neighbors in adjacent zig-zags, Ja. To explain our data it is necessary to allow 
for second neighbor interactions, Jf 2 and Jfs, within zig-zags (see main text). Nearest neighbor antiferromagnetic coupling 
between planes, Jc, is not shown, (b) The ZP model, indicating exchange interactions between rigidly coupled S = 7/2 units in 
the a&-plane. The ZP units are indicated by dashed ovals. Parallel units have AFM coupling, Ja, and perpendicular units have 
ferromagnetic couplings, either Jfu or Jfd. (c) The dimer model, as proposed by Johnstone et a^, in which spins are not 
rigidly coupled as in the ZP case, but rather have a stronger ferromagnetic intra-dimer coupling, Jfs, and weaker inter-dimer 
couplings, FM Jfw, and AFM Jai and Ja 2. The orbitals are shown to have a mixed character of a superposition of dy2_^2 
id^2_^2) and d^^2_^2 {d^y2_^2), following the conventions used in model calculations for extensions to the DDEX mode l^^’^° 
discussed in the text. As originally proposed the model was purely phenomenological without reference to an electronic model, 
and the orbitals should be considered purely schematic in this panel, (d) An alternative more general dimer model (see text) 
in which we distinguish between the Mn sites with different proportions of dy2_^2 (d^2_^2) and d^^2_^2 {d^y2_^2) character, 
indicated here by purple and light blue coloring respectively. The orbital character together with Goodenough-Kanamori rules 
lead to FM Jfw 2, and AFM Ja, different to the model shown in (c). 


to distinguish between the various scenarios that arise 
within DDEX models from an investigation of the spin 
waves. The measurement of the spin wave excitation 
spectra using neutron scattering directly yields informa¬ 
tion about the strength of the magnetic interactions in 
a material. In addition to the strength of interactions 


defining the energy of features in the dispersion relations 
such as overall bandwidth and sizes of gaps, the directions 
and distances of the interactions define the periodicity in 
wavevector of those features. These periodicities can be 
used to distinguish between the Goodenough model, the 
case of strongly bound dimers (which we refer to as the 
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’ZP model’) and weakly dimerized models (’dimer mod¬ 
els’), all of which exhibit the same magnetic structure 
but arise due to different sets of exchange interactions. 

As discussed briefly above, in two earlier studies of 
the spin waves in CE-type half-doped manganites^i*^, 
single-layered Lai/2Sr3/2Mn04 and pseudo-cubic 
Ndo.sSro.sMnOa, the lower part of the spin wave disper¬ 
sion was measured and found to be well-explained by the 
Goodenough model including second neighbor exchange 
Jf 2 along the FM zig-zags (see Fig. HJa)). The limiting 
ZP model of strongly bound dimers, that is, when the 
intra-dimer exchange is overwhelmingly stronger than 
any other interactions in the system, was dismissed. The 
authors argued that of the nearest-neighbor inter-dimer 
exchange interactions in the ZP model, only that along 
the AF stacking direction results in significant dispersion 
because the other two (labeled by Jfc/ within a chain 
and Jfd between chains in Fig. [I{b)) are frustrated, 
resulting in weak or no dispersion perpendicular to the 
AF stacking direction. This may be the case for the non- 
collinear structure proposed in ref. [l^ for bond-centered 
charge ordering on the basis of symmetry, in which the 
pairs of spins in adjacent ferromagnetic dimers along 
the zig-zag structure are orthogonal. However, for the 
case of the known CE-type magnetic structure^d^d^, 
there is no such a-priori reason to be confident that 
ipu and ipD are similar. The importance of measuring 
the full set of dispersion relations is illustrated in Fig. 
m Considering still the limiting case of strongly bound 
dimers. Fig. shows theoretical spin wave dispersion 
relations for the Goodenough model (dashed line) and 
ZP model (solid line)dd At low energies the dispersion 
relations near the zone center are qualitatively very 
similar. The main difference, however, is the absence 
of the higher branches for the ZP model, because the 
binding of spins in pairs halves the number of spin wave 
modes (see Section |TVl Discussion). 

In light of the discussion above of the origin of both 
the intra and inter dimer magnetic exchange along the 
zig-zags originating from double exchange in the DDEX 
description, a more realistic dimer model needs to allow 
for weak dimerization, that is, the case when the two 
exchange parameters are similar. In a recent time-of- 
flight (ToF) neutron spectroscopy study^ of the bilayer 
manganite Pr(Cao.9Sro.i)2Mn207 (PCSMO) the higher 
branches of the dispersion which are present in the Good- 
enough model but not in the ZP model were successfully 
measured. This enabled the authors to show that the 
Goodenough model clearly provides the better descrip¬ 
tion of the ground state in this material. Furthermore, 
the observed dispersion required both second neighbor 
interactions Jf 2 and Jfs along the FM zig-zags, which 
was ascribed to indicating some itinerant electron char¬ 
acter. The authors also considered an alternative weakly 
dimerized model, shown in Fig. [ijc). The model allows 
for alternating stronger and weaker nearest neighbor FM 
exchange along the zig-zags on a purely phenomenologi¬ 
cal basis. In PCSMO the differences between this dimer 



FIG. 2. Calculated dispersion relations for the Goodenough 
model (dashed red line) and ZP model (solid blue line). The 
parameters for the Goodenough model are those obtained as 
the global best fit in our analysis (see main text). The param¬ 
eters for the ZP model are J fu = 1.6 meV, Jfd = —0.11 meV, 
Ja = 0.11 meV and the inter-planar coupling is 1.05 meV. 
These parameters were chosen to show that there exists a set 
of parameters for the ZP model such that at low energies it 
is practically indistinguishable from the Goodenough model, 
but that the dispersion at higher energies serves as an effective 
discriminant. The x-axis labels denote high-symmetry posi¬ 
tions in the first Brillouin zone - G = (0,0,0), X = (1,0,0), 
y = (0,^0), Z= (0,0,1), M = (1,1,0), and L = (1,1,1). 
Note that we have chosen to index reciprocal space according 
to space group Pbnm, described in the main text. 


model and the Goodenough model could be resolved in 
favor of the latter by inspecting maps of the intensity 
^(Q, w) at a particular key energy, namely the top of the 
lower band of excitations. 

In this paper we present an extensive investigation of 
the spin excitation spectrum of Pro.sCao.sMnOa. This 
material was chosen since it is perhaps the canonical ex¬ 
ample of a CE-type half-doped manganite. Furthermore, 
the whole debate surrounding whether the Goodenough 
or ZP model is an appropriate description of the interac¬ 
tions in the ground state arose because of the pioneering 
work of Daoud-Aladine et on PCMO close to half¬ 
doping. We compare our time-of-flight inelastic neutron 
scattering measurements qualitatively and quantitatively 
to each of the Goodenough, ZP and generalised weak 
dimer model shown in Fig. [TJd). We consider in detail 
the dispersion relations and dynamical structure factor, 
S'(Q, w), elucidating the specific features in the data that 
allow us to distinguish between the models. Despite the 
complication of six-fold twinning we are able to show 
that the Goodenough model provides the best descrip¬ 
tion of our data. It is our intention that this detailed 
description of how inelastic neutron scattering data can 
be analyzed to distinguish between the different scenarios 
that can arise in DDEX models will serve as a reference 
for future work on manganites away from half-doping, 
for which the calculations indicate a greater instability 
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towards the formation of dimer-like ground states. 


II. EXPERIMENTAL METHODS 

The ToF neutron scattering measurements were per¬ 
formed on the ARCS chopper spectrometer at the Spalla¬ 
tion Neutron Source, Oak Ridge National Laboratorj*^. 
Data were collected with incident neutron energies (Ei) of 
35, 70, and 140 meV, with the instrument’s correspond¬ 
ing Fermi chopper frequencies chosen in order to give en¬ 
ergy resolution Ri 0.05Ei (full-width half-maximum) at 
the elastic line. The sample used for these measurements 
comprised a single crystal of mass 1.57 g, grown as de¬ 
scribed previously^. It was mounted in an aluminium 
can containing helium exchange gas. The can was then 
mounted on the cold finger of a closed-cycle refrigera¬ 
tor (CCR) and cooled to 5.5 K, at which temperature all 
measurements were performed. The sample was mounted 
with the a- and c-axes horizontal. For the data collected 
with Ei = 35meV and Ei = 140 meV the sample was 
held in a single orientation, with the c-axis of the sample 
parallel to the incident neutron beam. The data in these 
configurations were collected for 4 hours and 29 hours 
respectively. We hereafter refer to measurements taken 
in this configuration as ‘single-shot’ datasets. 

For the data collected with = 70 meV, we performed 
a single-shot measurement with the c-axis of the sample 
parallel to the incident neutron beam, in which config¬ 
uration we collected data for 14 hours. In addition we 
performed a measurement in which the sample c-axis was 
rotated about the vertical axis from parallel to perpen¬ 
dicular to the incident neutron beam, in 0.5° steps, with 
spectra recorded for ^ 1 hour at each orientation. The 
spectra at each orientation were then combined into a 
single large dataset using the Horace software^, allow¬ 
ing us to make a broad survey of the cross-section in 
the full 4-dimensional wavevector-energy space. We here¬ 
after refer to measurements taken using this method as 
‘multi-angle’ datasets. Such a broad survey is crucial for 
a complete understanding of this material since, unlike 
bilayer PCSMO for example, the magnetic interactions 
have appreciable strength in all three dimensions. In a 
ToF measurement of this nature energy transfer is cou¬ 
pled to the component of the wavevector Q parallel to 
the incident neutron wavevector kj. Only by performing 
measurements in multiple orientations (or alternatively 
with multiple incident energies) can one decouple Q and 
energy^. To illustrate this, we show in Fig. [3]the ARCS 
spectrometer’s detector coverage of the ac-plane of recip¬ 
rocal space for a single-shot run (c-axis parallel to the 
incident neutron beam) and for a 90° angle scan with 
Ei = 70meV. 

We have chosen to index reciprocal space according 
to space group Pbnm^. In this convention, the zig-zag 
chains propagate along the (1,0, 0)-direction, as shown 
in Fig. [Ha). This crystal structure is related to the 
cubic perovskite system, which has lattice parameter 
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FIG. 3. (Color online) (a) The reciprocal space coverage of 
the detectors in the (equatorial) {H, 0, 0)/(0, 0, L)-plane, with 
Ei — 70meV, zero energy transfer, and for a multi-angle 
dataset with the sample c-axis rotated stepwise from parallel 
to perpendicular to the incident beam, (b) Reciprocal space 
coverage for a single-shot dataset with the orientation fixed 
to be the c-axis parallel to the incident neutron beam. 

Op = 3.85A, by a ~ 6 ~ V^op, c ~ 2ap, thus the mag¬ 
netic unit cell has size 2-\/2ap x 2^/2ap x 2ap. We there¬ 
fore expect magnetic Bragg peaks, in the Pbnm conven¬ 
tion that we use, with wavevectors ki = (0,1/2, L) and 
k 2 = (1/2,1/2, A), and structural peaks with wavevec¬ 
tor q = (1/2,0, A'), where A is an odd integer and L' is 
any integer, arising as a result of a Jahn-Teller distortion 
of the Mn-0 octahedra and orbital ordering^. In PCMO 
we expect sixfold twinning of the crystal that is described 
by a systematic permutation of the principal axes of the 
pseudocubic lattice. In order to demonstrate the effect 
of the twinning we show in Fig. |T]the various peaks we 
expect to observe for a series of planes with A = 0, 0.5 
and 1 respectively. In the analysis of our INS data we 
will account for the twinning in the modeling explained in 
SectionHYl In addition, we expect to observe crystal elec¬ 
tric field transitions originating from the Pr ions. Inelas¬ 
tic neutron scattering experiments^ on Pro.sSro.sMnOg 
show that there are transitions with energies 12.8 meV 
and 15.35 meV. We did not observe strong transitions 
(compared to the other sources of magnetic scattering) 
at these energies, however. 


III. RESULTS 

Figure|5]shows an overview of the dispersion from mea¬ 
surements with all three incident energies used. The 
essential features shown here are that the dispersion is 
highly structured over the full energy range, but most 
importantly there are two bands of dispersive excitations 
with an energy gap between them. Qualitatively, we 
estimate the two bands to cover energies in the range 
0 < E <37 meV and 42 < A < 75 meV. The data shown 
in panels (a) - (c) are taken from single-shot datasets, 
so the component of Q parallel to the incident neutron 
beam, (0, 0, A), varies as a function of energy. This makes 
qualitative statements difficult to make, other than the 
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FIG. 4. (Color online) The magnetic Bragg peaks ki = 
(0,1/2, L) and k 2 = (1/2,1/2, L), and structural superlattice 
peak q = {1/2,0, L'), indicated by red circles, blue squares 
and green triangles respectively. Here we have taken account 
of the six-fold structural twinning of the PCMO crystal. No¬ 
tice that the k 2 peaks do not overlap with the structural 
superlattice peaks, whereas the ki peaks do. 
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obvious presence of two bands of spin waves, but the 
value of L is known for each energy so these data can 
be used collectively to constrain the parameters of differ¬ 
ent spin wave models. In panel (d) we show data from 
the multi-angle dataset, so the value of L is the same for 
all of the data presented. This serves to illustrate that 
the gap observed between the lower and upper bands of 
the spin wave dispersion is not an artefact of the way 
in which a gapless dispersion might project on to the 
curved (Q, E) hypersurface of measured reciprocal space 
coordinates from a single-shot run. 

To follow the evolution of the spin waves with energy, 
we show in Fig. |6] a series of intensity maps at fixed 
energy from single-shot runs, together with simulations 
using the Goodenough model with the global best fit pa¬ 
rameters obtained in the next Section. Because the value 
of L changes with incident energy and energy transfer for 
a given Q coordinate in these plots, we are mapping the 
L dispersion as well as the in-plane dispersion. The first 
four panels ((a) to (d)) therefore show a spin wave which 
disperses from (1,1), with the fourth panel, (d), show¬ 
ing the top of the lower band of the spin waves. Panels 
(e) and (f) show the upper band of spin waves, with the 
dispersion converging on (1, 0) as the band maximum at 
~ 75 meV is approached. This part of the spectrum can 
be understood, with reference in addition to the period¬ 
icity evident in Fig. EKc), as arising from dominant FM 
exchange between nearest neighbor Mn ions along the 
zig-zags. 

We incidentally note that the gross features described 
above are not dependent on a consideration of the struc¬ 
tural twinning that is present in this material. That is 
to say, the presence of a gap must be a feature of all the 
twins present in our sample, and is independent of their 
relative populations. Similarly, the periodicity of the ob¬ 
served signal in momentum space can only arise from 
the periodicity of the spin chains. This is as opposed 
to a situation where contributions from different twins 
that might each have different periodicity that taken to¬ 
gether combine to give an overall signal with periodicity 



FIG. 5. (Color online). Neutron scattering intensity maps of 
the dispersion along the {H, 0)-direction, which is parallel to 
the zig-zag chains. The maps are taken from the single-shot 
datasets with (a) Ei = 35meV, (b) Ei = 70meV and (c) 
Ei = 140 meV. The data have been smoothed by convolution 
with a Gaussian with FWHM equal to two bin sizes along 
each dimension for clarity. The map in panel (d) is taken 
from the multi-angle dataset with Ei = 70meV, hence why 
the value of L can be explicitly given. The white dashed 
rectangles in panels (b)-(d) indicate the position of the gap 
between upper and lower branches of the excitations, which is 
consistent with the Goodenough model. Note that the data 
in all four panels have been multiplied by f{E) = ^ i 

where Eq = ks x lOK. f{E) is therefore constant when E <C 
Eo and f{E) oc E when E ^ Eq. This rescaling is to display 
the full spectrum more clearly on a single intensity scale. For 
panel (a) the signal was integrated over ±0.05 r.l.u. in the 
{0,K) direction, whereas for panels (b)-(d) the integration 
was ±0.1 r.l.u. 


not seen in any single twin. 


IV. DISCUSSION 

The DDEX mode l^^d^ , which is a more general ver¬ 
sion of double exchange, would be an appropriate start¬ 
ing point for analyzing our data. However, calculation 
of the excitation spectra within such a framework is, 
though possible, far from straightforwar d^^i^° . In the 
absence of such calculations, we therefore use an effec¬ 
tive Heisenberg Hamiltonian to model our data, just as 
was done for the half-doped manganites mentioned in 
the introductio n^^d^d^ .The use of an effective Heisen¬ 
berg Hamiltonian is a well-established procedure to an¬ 
alyze the spin wave spectra in magnetically ordered sys¬ 
tems, and not just those with purely integer or half¬ 
integer spins, nor even solely insulators. The exchange 
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FIG. 6. (Color online). Panels (a) - (f) show intensity 
maps taken with energies in the range 9<i?<12, 15<i?<20, 
20<i?<25, 35<i?<40, 50<i?<60, and 70<i?<80meV, and in¬ 
cident energies of 35, 70, 70, 70, 140, and 140 meV respec¬ 
tively. The data have been smoothed by convolution with 
a Gaussian with FWHM equal to two bin sizes along each 
dimension for clarity. Panels (g) - (1) show corresponding 
simulations of those data using the Goodenough model and 
the global best fit parameters (see main text). All of the data 
shown here are from single-shot datasets with the c-axis par¬ 
allel to the incident beam. We have made use of the symmetry 
of the crystal, and the symmetry and extent of the detector 
array by folding the data about the (1, 0, 0) and (0,1, 0) axes 
to effectively increase the overall count rate by a factor of four. 
The calculated S{(^,lj) were convoluted with Gaussians with 
FWHM matched to the spectrometer’s energy resolution at 
each energy and Ei . 


interactions in the model Hamiltonian that are required 
to reproduce the observed dispersion relations and in¬ 
tensities establish the strength, range and symmetries 
of the magnetic interactions. Examples where such 
an analysis is done include the parents of iron based 
superconductors^ii^ in which the signs and magnitudes 
of the nearest neighbor and further neighbor exchange 
constants elucidated the dominant second neighbor an¬ 
tiferromagnetic exchange. In the case of manganites ex¬ 


amples include the ferromagnetic metallic-like materials 
with X ~ 0.2 — 0.48 (e.g. refs. where softening of 

the dispersion towards the zone boundary and the doping 
dependence of intra-bilayer exchange could be qualita¬ 
tively explained by mapping the excitations of the double 
exchange model onto an equivalent ferromagnetic Heisen¬ 
berg Hamiltonian (valid at the quasi-classical level) with 
corrections beyond 1/S and incorporating the Cg orbital 
occupancies^i^i^. 

To analyze the data quantitatively we calculated the 
dispersion relations and dynamical structure factors for 
the Heisenberg Hamiltonian TL = — Jij®i ’ ~ 

D(Sf)^ using linear spin wave theory for each of the 
Goodenough, ZP and dimer models illustrated in Fig. 
[T](a-d). In the present analysis appropriate choice of ex¬ 
change interactions yields the known magnetic ground 
state for all of the models. The procedure for calculat¬ 
ing the dispersion relations for the first two cases is ex¬ 
plained in refs, andls^. in which textbook linear spin 
wave theory^ is applied to this particularly complex ex¬ 
ample. An extended explanation, including calculation 
of the structure factors, and will be the subject of a fu¬ 
ture publication. For the generalized dimer model the 
linear spin wave theory dispersion relations and dynam¬ 
ical structure factors were calculated using the McPhase 
package^. We conhrmed that in the special case when 
this model reduces to the Goodenough model that it re¬ 
produces the results of ref. [5^ In the first part of the 
analysis procedure for each model, we obtained robust 
initial estimates for exchange constants from examination 
of the periodicities and energies of key spectral features 
detailed below. We then fitted the calculated 5'(Q,u;), 
averaged with equal weight from each twin, to the data 
by least-squares refining the exchange constants Jy. The 
term D is a single-ion anisotropy term to constrain the 
spins to the a&-plane. This term was included with a 
fixed value of 0.07 meV in all of the models discussed in 
our analysis, to be consistent with the literatur o^^d^H^ ^ 
Its small value means that it has negligible effect on the 
cross-section at the energies that we probed in our ex¬ 
periment and hence does not affect any of the results of 
our fits or our conclusions. The data that we fitted were 
all either one-dimensional cuts or two-dimensional slices 
taken at constant energy, so that the energy-dependent 
resolution of the spectrometer could be accounted for by 
using a fixed Gaussian broadening of the calculated in¬ 
tensity, the width of which was chosen to match the reso¬ 
lution at each energy. Twenty eight one-dimensional cuts 
and twelve two-dimensional slices were used in the fitting 
process, with the fit being performed simultaneously on 
all of them to determine a single set of best-fit parame¬ 
ters. The robustness of the fit parameters obtained was 
checked by using numerous different sets of starting val¬ 
ues, to confirm that the fit always converged on the same 
result. 

At this point we reiterate the distinction between the 
different models we will consider. In all cases they have 
the same magnetic ground state, as shown in Fig. [TJ 







































which is commonly referred to as the CE magnetic struc¬ 
ture. What distinguishes what we refer to as the ‘Good- 
enough model’ and the ‘ZP model’ and dimer models is 
thus the nature of the electronic environments around 
the Mn ions, which give rise to different sets of ex¬ 
change interactions between them. Our INS measure¬ 
ments, and analysis with the effective Heisenberg model, 
therefore provide insight into the nature of this aspect 
of the ground state, rather than the magnetic structure 
itself, which is well establishedi^. 


A. Goodenough model 

We start with the Goodenough model. We use this 
term to apply to a model which has the CE magnetic 
structure and a set of exchange parameters as shown in 
Fig. [ija). Although we consider during our analysis the 
possibility of charge disproportionation, we do not a pri¬ 
ori assume that it is present. The orbital order is implic- 
ity assumed to be that shown in Fig. [Ija), since this then 
provides a natural explanation of how the signs of the 
nearest neighbor exchange interactions we consider come 
about. We also allow for the possibility of next-nearest- 
neighbor interactions along the EM zig-zags, as has been 
required in previous studies of similar system o^^did^ and 
is taken to indicate a degree of electronic itinerancy along 
the zig-zag chains. 

To simplify the analysis we initially assumed that the 
charge disproportionation between notional Mn^“*" and 
Mn^+ sites is zero. This meant that all spins were set to 
an average value of S' = 7/4. It is well known that the 
charge disproportionation is substantially smaller than 
the ideal value^ d^d^ and close to Mn^-^+ and Mn^ ®+. 
By choosing to use equal spins on all sites, we are able to 
draw conclusions more straightforwardly about the mag¬ 
nitude and sign of exchange interactions from qualitative 
analysis of the periodicity and extent of the dispersion. 
Later in our analysis we will revisit this assumption and 
will show that the effect of introducing a magnetic mo¬ 
ment disproportionation equal to the experimentally de¬ 
termined valuei^ does not alter the principal conclusions. 

As has already been noted, the dominant exchange is 
the nearest-neighbor interaction ipi linking adjacent Mn 
ions along the zig-zag chains. Examination of Fig. [Sjc) 
shows that the maximum of the dispersion, which has the 
periodicity of a ferromagnet, is ^ 75 meV, and thus the 
strongest constraint on the size of J pi comes from this 
part of the dispersion. We determined from our global 
best fit that Jfi = 9.01 ± 0.02 meV. This value is compa¬ 
rable to that measured in PCSMO (11.39 ± 0.05 meV)^, 
Lao. 5 Sri. 5 Mn 04 (9.98 meV)^, the AFM CO phases of 
Ndo.sSi'o.sMnOs (6.76 meV)^, and a little larger than 
the typical values for the EM double exchange in metal¬ 
lic EM manganites22i^. In other words, it is consistent 
with its origin being due to double exchange. 

At much lower energies (^ lOmeV) there is a branch 
of the dispersion, shown in Fig. [71 with double this peri¬ 


odicity which arises from the antiferromagnetic coupling 
between zig-zag chains, J^. Variations in this parameter 
have negligibly small effects on features of the dispersion 
above ~ 10 meV, so in practice it can be determined from 
fits to the Ei = 35meV dataset alone, where the back¬ 
ground arising from the tails of the elastic line at these 
low energies is minimized. The global best fit yielded a 
value of Ja = —0.45 ± 0.01 meV. 



-0.5 0 0.5 

(0.5,K,1) [r.I.u.] 


FIG. 7. (Color online). Low energy dispersion, the extent 
in energy of which allows a reasonable estimate of J a to be 
obtained before the global fitting procedure. 

The existence of a gap in the dispersion at (1/2,0) 
(Fig. |5l(d)) immediately implies a halving of the period¬ 
icity in real space along the zig-zags. In the framework 
of an effective Heisenberg Hamiltonian this implies one 
or both of the second neighbor interactions Jf 2 and 
is non-zero. The position in energy and size of the gap 
between the lower and upper spin wave bands, shown 
schematically in Fig. |2l define two energies which deter¬ 
mine the values of both the next-nearest-neighbor terms 
Jf 2 and Jfs- In order to reproduce the measured values 
Jf 2 and Jf 3 must both be non-zero, and we find that 
in fact they must have opposite signs, as was similarly 
concluded for PCSM022. The size and position of the 
gap alone can be reproduced identically by exchanging 
the signs of Jf 2 and Jf 3 , that is, with Jf 2 < 0 <Jf 3 
or Jf 2 > 0 >Jf 3 - However the details of the structure 
factor around the gap can distinguish between these two 
scenarios, specifically spectral weight is swapped between 
two relatively flat branches of the dispersion that lie at 
the top of the lower set of spin wave bands and the bot¬ 
tom of the upper set of spin wave branches. Also, the 
periodicity of peaks in the structure factor at the top of 
the lower set of spin wave branches is doubled along a 
key symmetry direction for the case of Jf 3 > 0 >Jf 2 , at 
variance with the data. This is illustrated in Fig. |8l For 
the global best fit we find that Jf 2 = 0.93 ± 0.02 meV 
and Jf 3 = —0.97 ± 0.02 meV. Note that Jf 3 is therefore 
antiferromagnetic. 

A particular difference between PCMO and PGSMO 
is that in the former the coupling between planes is ex- 
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FIG. 8. (Color online). Cut across part of the top of the lower 
band of excitations. The red (solid) line shows the global best 
fit to the data in this region, whereas the blue (dashed) line 
shows the structure factor if the signs of both Jf 2 and Jf 3 
are reversed. Such a reversal preserves the magnitude of the 
gap between the lower and upper bands of excitations, and 
also the dispersion relation, but changes the structure factor 
as shown. 


pected to be significant, whereas in the latter it is rather 
small. Indeed, previous low energy INS measurements 
of Ndo.sSro.sMnOa (ref. EH), which one would expect to 
be similar to PCMO, found that the inter-plane coupling 
was about a factor of four larger than the coupling be¬ 
tween zig-zag chains. The complex twinning, together 
with an intermediate strength coupling between planes, 
means there is not a single feature which can be used 
to determine Jc easily. The strength and position of the 
signal around (±0.5,0) at ~ 22meV in Fig. [IKa), with 
Ei = 35meV, exhibits qualitative dependence on Jc, but 
overall the best method for refining Jc is to use all of 
the data when performing the final global fit. We found 
Jc = —2.09 ± 0.03 meV, which like in Ndo.sSro.sMnOa is 
about a factor four larger than the in-plane antiferrog- 
magnetic coupling between zig-zag chains. 

Let us now return to the question of charge / spin 
disproportionation between the different Mn sites. The 
effect of this being non-zero would be to have different 
values of the spin on alternate Mn sites along the zig¬ 
zags, which would halve the real-space periodicity and 
thereby open up a gap to create lower and upper disper¬ 
sion branches, even with Jf 2 =Jj ’3 = 0. 

To account for different moments on the Mn sites we 
can set S = 1.75 ± SS on the notional Mn^+ sites and 
S = 1.75 — 5S on the notional Mn^+ sites, instead of 
the uniform values of 5" = 1.75 used up to now. In 
fact, the dispersion relations throughout the Brillouin 
zone for our global best fit parameters for the Goode- 
nough model can be almost identically reproduced with 
Jf 2 =Jf 3 = 0, but this requires SS = 0.42. This corre¬ 
sponds to a ratio of moments on the Mn^’*' / Mn^’*' sites 
of 1.63, significantly at odds with the measured ratio^^ 


of 3.18/XB/2.75/iB = 1.16. Alternatively, if the differ¬ 
ence in the spins is considered to arise from charge dis¬ 
proportionation in a purely ionic model (in which Mn^’*' 
has S = 2 and Mn'^^ has S = 3/2) this corresponds 
to Mn^-^+ and Mn^-^+. In contrast, bond valence sums 
for PCMO from joint refinement of neutron and x-ray 
diffraction dataM yield Mn^ '^®+ and Mn^ ®^+, consis¬ 
tent with typical disproportionation of ± O.le in half- 
doped manganites^ii^rii. The gap therefore cannot be 
accounted for by the experimental moment or charge 
disproportionation. The contribution that it can make 
can be determined as follows. If we fix the ratio of the 
spins to that of the measured magnetic moments, then 
this changes the values of J_f 2 and J_f 3 to 0.70 meV and 
—0.64meV c.f. 0.93 meV and —0.97meV respectively, a 
reduction in the values of about 30%. Alternatively, fix¬ 
ing the ratio of the spins by the values corresponding to 
typical charge disproportionation of ±0.1e changes 3f2 
and Jfs to 0.84 meV and —0.84meV, a reduction com¬ 
pared to the values assuming no disproportionation of 
less than 14%. To summarize, the difference between the 
Mn moments or charge disproportionation on the two 
sites can therefore only account for a small fraction of 
the measured gap. The position in energy and the size of 
the gap therefore robustly requires non-zero J f 2 and J ^3 
in the Goodenough model. We additionally note that ac¬ 
counting for the small charge / spin disproportionation 
as discussed above only results in a change of the values 
of 3f2 and J^s ; the other fitted exchange parameters, 
Jj’i, 3a and Jc, remain unaltered within the size of the 
errorbars. 

Our analysis has up to this point has also assumed 
an equal population of the six structural twins that are 
possible in this system. It is simplest to describe the 
twins with respect to the pseudo-cubic perovskite unit 
cell. Twin 1 is the unrotated coordinate system, and twin 
2 corresponds to a 90° rotation of this system about its 
pseudo-cubic c-axis. Twin 3 is a 90° rotation of twin 1 
about its a-axis, and then twin 4 is equivalent to twin 3 
rotated by 90° about the c-axis of twin 1. Twin 5 is a 
90° rotation of twin 1 about its 6-axis, and finally twin 
6 is a further 90° rotation of twin 5 about the c-axis of 
twin 1. 

We now illustrate the validity of the approach of as¬ 
suming an equal population of twins in two ways. First, 
in Fig. [5] we show simulations of some of the constant- 
energy slices that have already been shown in Fig. |6l 
Specifically, panels (a) and (b), and (c) and (d), should 
be respectively compared to panels (g) and (h) of Fig. 
m and panels (e) and (f) should be compared to panel 
(1) of Fig. [6l We have re-calculated the scattering in¬ 
tensity for the global best fit to the Goodenough model, 
decomposed into the contribution from each of the six 
twins, and then formed combinations of them to show 
that there is no permutation of fewer than all six twins 
that will allow us to reproduce the fits to our dataset 
across the entire energy range. For example, in Fig. [91(a) 
we can qualitatively reproduce the global best ht with 
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FIG. 9. (Color online). Decomposition of some of the fits pre¬ 
sented in Fig. [6]into contributions from specific sets of twins, 
(a) Contribution from twins 1-4, which should be compared 
to Fig. o a) and (g); (b) contribution from twins 5 and 6, 
i.e. those which are not included in panel (a); (c) contribu¬ 
tion from twins 1, 2, 5, and 6, which should be compared to 
Fig. [6Kb) and (h); (d) contribution from twins 3 and 4, i.e. 
those not included in panel (c); (e) contribution from twins 
3-6, which should be compared to Fig. (SJf) and (1); (f) con¬ 
tribution from twins 1 and 2, i.e. those not included in panel 
(e). Our convention for labeling the twins is described in the 
main text. 


contributions from the twins we label 1 — 4, whereas in 
Fig. iKc) we require contributions from twins 1, 2, 5 
and 6, and in Fig. j^Ke) all of the features in the data 
can be accounted for with just twins 3 — 6. Second, we 
repeated the entire fitting procedure with additional fit 
parameters of the relative contributions from each twin. 
We found that, irrespective of the initial values of the fit 
parameters chosen, within the error on the fit parameters 
the populations of all six twin were essentially equal. 


B. Zener polaron model 

Let us now turn to the ZP model (that is, of strongly 
bound dimers) for the sake of completeness, before con¬ 
sidering the more realistic dimer model in Sec. IIV Cl As 
shown in Fig. ^ a key feature that distinguishes the 
Goodenough and ZP models is the presence or absence 
of a second band of excitations. By far the dominant ex¬ 
change constant in the ZP model is expected to be that 
within the dimers. Accordingly, we expect one set of spin 
wave branches associated with coupling of the rigid ZP 
dimers that will have a band width determined by the 
characteristic inter-dimer exchange energies. These will 
be separated from a second set of branches associated 
with excitation of the pair of spins within the ZP by an 


energy gap determined by the much stronger intra-dimer 
exchange. That is, we expect a gap that could be sev¬ 
eral multiples of the bandwidth of the lower energy set 
of spin wave branches, which in the limiting case of the 
ratio of intra-dimer to inter-dimer exchange tending to 
infinity eliminates half the spin wave branches. In con¬ 
trast, the Goodenough model permits a small (or even 
zero) gap between two sets of branches depending on the 
existence (or otherwise) of second neighbor interactions 
along the ferromagnetic zig-zag chains. As illustrated in 
Fig. m judicious choice of Ja, Jfu and Jfd in this lim¬ 
iting ZP model (the model shown in Fig. (Ub)) results 
in spin wave branches that closely follow the lower en¬ 
ergy half of the spin wave branches of the Goodenough 
model, but in general the reverse is not possible. This is 
because the limiting ZP model has only half the number 
of spin wave modes by virtue of the rigid coupling of spins 
within the dimers. Moreover, the periodicity of the most 
energetic branch in our data corresponds to the nearest- 
neighbor Mn-Mn distance projected along the direction 
of the zig-zags. The nearest-neighbor separation of ZP 
dimers is twice that, which would result in a doubling 
of the periodicity of the spin wave branch compared to 
what we observe. In summary, because we have been able 
to measure the higher energy branches, which are sepa¬ 
rated from the lower branches by a relatively small gap, 
we are able to discount a model of weakly interacting 
dimers without making assumptions about the relative 
magnitudes of Ja, Jfu and 3fd- In fact, this conclusion 
should apply to any case of interacting dimers, so long as 
the intra-dimer exchange is dominant over all other inter¬ 
actions (nearest neighbor and further neighbor) between 
a spin in one dimer and those in other dimers. 


C. Dimer model 

As discussed earlier, it is more realistic to consider 
a model of weak dimerization as an alternative to the 
weakly interacting strongly bound dimers of the ZP 
model. A model along these lines was discussed by John¬ 
stone et al in their work on PCSMO^, and is shown in 
Fig. HKc). The empirical rationale for this model was 
that the dominant exchange is the FM nearest-neighbor 
interaction along the zig-zags, but allowing for a small 
modulation between weaker and stronger exchange will 
create a gap at (1/2,0) by virtue of doubling the repeat 
distance along the zig-zags. In PCSMO, with a suitable 
choice of parameters, this model produces an essentially 
indistinguishable dispersion relation to the Goodenough 
model, in addition to the accepted GE magnetic struc¬ 
ture. However, the two models can be distinguished by 
differences in the structure factor at the top of the lower 
band of the dispersion. 

The weak dimer model of ref. [s^ was purely phe¬ 
nomenological without reference to an electronic model. 
However, it is in fact a special case of a more general 
weak dimer model, shown in Fig. (IKd), which we will use 
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FIG. 10. (Color online), (a) Intensity map of the dispersion 
along (—1.5, if, 1), taken from the Ei = 70meV multi-angle 
dataset. The energy-dependent background arising from the 
tail of the incoherent elastic line has been subtracted. The 
signal was integrated over ±0.1 r.I.u. in the perpendicular 
{H, 0, 0) and (0, 0, L) directions, (b) Simulation of the data 
in panel (a) using the global best fit parameters of the Good- 
enough model, (c) Simulation of the data in panel (a) us¬ 
ing the global best fit parameters for the weak dimer model. 
These parameters give an essentially identical dispersion to 
the Goodenough model, with the only difference being the rel¬ 
ative spectral weight at (—1.5, 0, l)-type positions compared 
to that at (—1.5, 0.5, l)-type positions at the top of the lower 
spin wave band. 


in the following discussion. The electronic justification 
for alternating stronger and weaker nearest-neighbor FM 
exchange J ps a-nd J fwi arising from double exchange 
along the zig-zags has already been described in detail 
in the introduction. For the time being we disallow the 
next nearest neighbor exchange Jf 2 and Jf 3 along the 
zig-zags, but will include them later in this section. We 
also set all spins to an average value of S' = 7/4, as 
we did for the Goodenough model. In order to antic¬ 
ipate the likely signs of the weaker superexchange ex¬ 
change interactions between zig-zags we need to con¬ 
sider the orbital character on the Mn sites. Model cal¬ 


culations for extensions to the DDEX Hamiltonia n^ ^ 
that include on-site Coulomb repulsion, MnOe tilting and 
Jahn-Teller distortions, consider the Mn orbitals to be a 
superposition of dy 2_^2 [dx' 2 -z 2 ) and {d^y 2 _x 2 ) 

character, where the x and y directions point along the 
two perpendicular Mn-Mn directions along the zig-zags. 
When the parameters of the model yield dimers there 
is a qualitatively closer resemblance to the latter on all 
sites. At the corner sites along the zig-zags the resem¬ 
blance is less pronounced and tends to equal character 
(i.e. no orbital ordering) as the dimerization vanishes. 
Using the Goodenough-Kanamori rule a^°i^° for the ex¬ 
pected superexchange between zig-zags, we might expect 
AFM coupling between nearest-neighbor Mn ions in par¬ 
allel dimers, and FM coupling between nearest-neighbor 
Mn ions in perpendicular dimers. We therefore denote 
the AFM interaction between nearest-neighbors in par¬ 
allel dimers with the exchange parameter J^, and the 
weak FM interaction between nearest-neighbors in per¬ 
pendicular dimers with the exchange parameter 3fw2- 
This model differs from that proposed by Johnstone et 
al, in which it was assumed 3a2 =3fw2 < 0. Note that 
although we use the Goodenough-Kanamori rules to an¬ 
ticipate the signs of Ja and 3fw2, we do not actually 
constrain the signs in our modelling. Consequently, the 
dimer model shown in Fig. [TJc) is a limiting case of the 
more general model shown in Fig. m- 

We htted the dimer model in the same manner as the 
Goodenough model. An estimate for the strong intra¬ 
dimer coupling 3FS was made from the highest-energy 
part of the dispersion. The size of the gap between 
the upper and lower spin wave branches is mostly de¬ 
termined by the difference between the intra- and inter¬ 
dimer couplings along the zig-zags, which enabled an es¬ 
timate for J FWi to be made. Estimates for the weaker 
couplings between the zig-zags, 3a and 3fw2: were de¬ 
termined from low-energy features of the data. The 
dimer model parameters were then refined using the same 
cuts as were used for the refinement of the Goodenough 
model parameters. The global best fit parameters were 
then J FS = 10.74 ± 0.03 meV, J fwi = 7.02 ± 0.04 meV, 
3fw2 = 1.26 ± 0.11 meV, = —2.47 ± 0.14 meV, and 
Jc = —2.19 ± 0.03 meV. This produces an almost identi¬ 
cal dispersion relation to that obtained from the global 
best fit parameters to the Goodenough model. 

The small size of the gap between the spin wave 
branches compared to the overall bandwidth is reflected 
by rather weak dimerization, quantihed by the difference 
between J fs a-nd J fwi as a fraction of their average, be¬ 
ing only 42%. It is therefore not surprising that both J fs 
and J FWI are quite close in value to J fi in the Good- 
enough model, since they are largely determined by the 
maximum energy of the dispersion. As an interesting 
aside, the magnitude of the AEM coupling between the 
zig-zag chains in the best fit to the dimer model is a factor 
~ 5.5 larger than in the Goodenough model. This can 
be understood by noticing that is used for the two 
exchange pathways in the Goodenough model that are 
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described by and Jfvc 2 in the dimer model. In order 
for the same magnetic structure to be stable in the dimer 
model, |Jyi| has to be larger than Jfvk 2 , since there are 
two AFM contributions (J^) to the total energy in the 
Goodenough model, balanced against one AFM contribu¬ 
tion (Ja) and one FM contribution (3fw2) in the dimer 
model. These last two terms effectively trade against one 
another to give an overall antiferromagnetic interaction 
between zig-zags and close to identical dispersion to that 
observed in the Goodenough model. The interaction be¬ 
tween the planes that contain the zig-zag chains, Jc, is 
almost the same for the dimer and Goodenough models, 
since this term describes the same physics in both cases. 



FIG. 11. (Color online). Cut throngh the data shown in Fig. 
[inia), integrating the signal over 29 < E < 34meV. The 
data are indicated by black circles. The simnlations using the 
Goodenough and dimer models are indicated by dashed (red) 
and solid (blue) lines respectively. 

One might be tempted to conclude that because the 
dispersions obtained from global best fits to both the 
Goodenough and dimer models are essentially indistin¬ 
guishable, the models are equivalent. In fact this is not 
the case. For most of the energy scale of the excitations 
the two models have very similar structure factors. How¬ 
ever we find that near the top of the lower spin wave band 
the structure factors are different, and are readily distin¬ 
guished on examination of the data along a well-chosen 
symmetry direction. This difference in structure factor 
between the data and Goodenough model, and the dimer 
model, is illustrated in Fig. 1101 The agreement between 
the data and Goodenough models (panels (a) and (b)) is 
good, but there is enhanced intensity around (0.5, K, 1) 
for integer K arising from a larger structure factor from 
the dispersion which converges on these points in the 
dimer model (panel (c)) which does not agree with the 
data. To quantify this difference further we show in Fig. 
Illl a one-dimensional cut across the top of the lower band 
of the dispersion. The dashed line is the Goodenough 
model global best fit, which clearly captures all of the es¬ 
sential features of the data, whereas the solid line is the 


dimer model global best fit, that fails to do so. The differ¬ 
ences are clearly statistically significant, since the differ¬ 
ence between the dimer model and the data at (0.5, AT, 1) 
for integer K is much larger than the size of the errorbars 
on the data points. 

The essential physical difference between the dimer 
and Goodenough models is that the former allows for the 
gap at (1/2,0) and ~ 40meV to be accounted for by al¬ 
lowing alternating nearest-neighbor FM exchange along 
the zig-zags, whereas the latter accounts for it using fur¬ 
ther neighbor terms along the zig-zags. A question which 
naturally arises is what is the result if both weak dimer¬ 
ization and further neighbor interactions are permitted. 
The absence of peaks in the structure factor at (0.5, AT, 1) 
for integer AT shown in Fig. [TT] strongly supports the pure 
Goodenough model without any element of dimerization, 
as the peaks appear to be a feature solely of the dimeriza¬ 
tion. We can quantify this assertion by fitting the data 
in Fig. 111! with a generalization of the dimer model that 
adds the same further neighbor terms that were used in 
the Goodenough model, Jf 2 and Jfs- With the same ex¬ 
change interactions between the FM zig-zags, the disper¬ 
sion and structure factor of the Goodenough model are 
reproduced exactly in such a dimer model with the alter¬ 
nating nearest-neighbor FM exchange terms set to the 
same value i.e. zero dimerization, and the further neigh¬ 
bor terms equal to the Goodenough model values. For 
the fit, we fixed the overall bandwidth of the spin waves 
(which fixed the average of nearest neighbour FM ex¬ 
change, (Jfs+Jfwi)/2)), the energies of the centre and 
size of the gap (which further constrains Jf2 = — Jfs), 
and the inter-chain exchange terms, to the values from 
the global best fit dimer model, and then allowed the 
dimerization (the difference between J ps and J fwi as 
a fraction of their average) and Jf2(= — Jfs) to vary. 
The result is Jf 2 = — Jf 3 = 0.84 ± 0.09 meV, very 
close to the values for the global best ht to the Good- 
enough model in Sec. IIV Al (Jf 2 = 0.93 ± 0.02 meV and 
Jf 3 = —0.97 ± 0.02 meV), and the dimerization is re¬ 
duced from 42% in the dimer model without the further 
neighbor terms to 6.5 ± 6.6%. We therefore conclude 
that adding further neighbor exchange interactions to the 
dimer model permits the best fit to revert essentially to 
the Goodenough model. 


V. CONCLUSION 

We have studied the spin waves in PGMO using time- 
of-flight neutron spectroscopy and found that the upper 
and lower branches of the dispersion are gapped. The 
gapped spin wave dispersion is consistent with the Good- 
enough model with strong nearest-neighbor ferromag¬ 
netic interactions along the zig-zag chains and weak anti¬ 
ferromagnetic interactions between them, provided there 
exist small next-nearest-neighbor couplings Jf 2 and Jf 3 
along the zig-zags. Gross features of the spin wave dis¬ 
persion allow us to rule out the Zener polaron picture 
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in the limit of strongly bound dimers as a description 
of the ground state of PCMO. An alternative model in 
which pairs of Mn ions are weakly dimerized, is shown to 
provide an equally good description of the spin wave dis¬ 
persion relation as the Goodenough model, but there are 
significant differences between the spin wave intensities 
in this weak dimer model and the data. These differences 
do not arise in the Goodenough model, so our data indi¬ 
cate that the Goodenough model (an undimerized model) 
with significant second neighbor interactions within the 
FM zig-zags provides the most satisfactory description of 
the magnetic ground state of PGMO. The existence, sizes 
and signs of J f 2 and Jf 3 are rigorously set by the gap in 
the spin wave dispersion. Any charge / moment dispro¬ 
portionation between the Mn sites along the zig-zags is 
only able to account for a small fraction of the gap. 

We found that the nearest neighbor FM interaction 
along the zig-zags is comparable or larger than the 
typical values found for the FM exchange in metallic 
manganites^i^, a result consistent with previous studies 
of other half-doped manganites^i^i^. The magnitudes 
of the second neighbor interactions along the zig-zags are 
about twice that of the antiferromagnetic exchange be¬ 
tween zig-zags. The significant further neighbor interac¬ 
tions are suggestive of delocalization of the Cg electrons 
along the zig-zags. These results are consistent with the 
theoretical analysis of the degenerate double exchange 
(DDEX) model for half-dopin g^'^d^ . In this model, a 
cooperative ordering of the charge and degenerate Cg 
orbital degrees of freedom takes place to accommodate 
the frustration between delocalization (which favors fer¬ 
romagnetism) and the antiferromagnetic superexchange 
between the core t 2 g moments. The result is FM zig¬ 
zag chains that are antiferromagnetically coupled via su¬ 
perexchange, but with the FM exchange along the zig¬ 
zags arising from the double exchange mechanism. This 
explains the fact that the nearest neighbor FM exchange 
we measure in PCMO is comparable in magnitude to that 


found in metallic manganites. In the DDEX model the 
double exchange overcompensates for the antiferromag¬ 
netic superexchange along the zig-zags, and the delocal¬ 
ization of the Cg electrons plausibly accounts for the sec¬ 
ond neighbor exchange along the zig-zags that exceeds 
the antiferromagnetic superexchange between the core 
t 2 g Mn moments. For x = 0.5 the charge modulation 
of < 0.1 e is site centered so that all Mn-Mn bonds along 
the zig-zag are identical. Extensions to the DDEX model 
that couple the electronic and lattice degrees of freedom 
allow for cooperative distortion or tilt of the MnOe units 
that can result in magnetic dimerization for a range of 
parameters in the model. However, our results indicate 
that any such coupling does not affect the magnetic ex¬ 
change interactions in PCMO, and that the magnetic de¬ 
grees of freedom are best described by the Goodenough 
model within the DDEX picture. 

Given that all of the toy models we have considered 
provide a way to examine different ground states that 
come out of the DDEX model and extensions to it, our 
analysis can serve as a template for how to understand 
future experiments on other manganites away from half 
doping, where the predictions from DDEX suggest that 
the CE / Goodenough picture is less likely to apply. 
Given that the balance between the different ground 
states in these models is rather delicate, having a toolkit 
available which can distinguish between them in real ma¬ 
terials with different doping will provide valuable input 
to further calculations. 
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